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Abstract 

We describe several techniques for the calculation of multi-loop integrals and their 
application to heavy quark current correlators. As new results, we present the four- 
loop correction to the second and third physical moment in the low-energy expan- 
sions of vector, axial-vector and scalar quark current correlators. Using a Ward 
identity, we obtain the third and fourth moment for the pseudo-scalar correlator. 
We briefly discuss the impact of these results on the determination of the charm 
quark mass and the strong coupling constant using lattice simulations for the cur- 
rent correlators and of the charm- and bottom-quark mass from experimental data 
for o~{e + e~ — > hadrons). 

Key words: Perturbative calculations, Quantum Chromodynamics, Dispersion 
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1 Introduction 



One of the phenomenologically most interesting applications of heavy quark 
current correlators is the determination of fundamental parameters of QCD 
via sum rules. Low-energy moments of the vector correlator can be com- 
pared to weighted integrals over the experimentally measured i?-ratio, R(s) = 
<j{e + e~ — > hadrons)/<r(e + e~ — > to determine values for the charm- 

and bottom-quark masses [T1l2ll3ll4"] . Recently "data" from lattice simulations 
have been used instead of experimental results and indeed this method has 
become a competitive way of extracting the charm-quark mass and the strong 
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coupling constant [5]. Although the correlators of all four (axial- vector, vec- 
tor, scalar and pseudo-scalar) currents can be used in the lattice method, the 
most accurate predictions presented in [5] were based on the pseudo-scalar 
correlator. 

In order to obtain precise values for the quark masses, it is mandatory to cal- 
culate higher order QCD corrections. To match the experimental precision this 
means calculations in four-loop approximation have to be performed. Improve- 
ments in both computer power and the techniques of multi-loop calculations 
have lead to significant progress in recent years. 

After the second moment at order a 3 of the vector current became available 
it was possible to construct a Pade approximant of the vacuum polarization 
function [6]. This also leads to a prediction for the higher moments. With the 
explicit calculation of the third moment this prediction can be checked and it 
becomes possible to further improve the Pade approximation [7|. 

The three-loop, i.e. 0(a 2 ), corrections were evaluated more than ten years ago 
in [8ll9lfT0] for vector, axial-vector, scalar and pseudo-scalar currents in the low 
energy expansion up to (q 2 ) 8 . This calculation has employed the reduction 
method proposed in [11] . Recently, a different approach based on the combi- 
nation of the Laporta algorithm [12] with differential equations [T3"|14lll5j was 
used to calculate these corrections for terms up to (g 2 ) 30 [TBfTT] . 

At four-loop order, for the vector correlator the first physical moment propor- 
tional to {q 2 ) 1 was obtained in [T8irT9] . The second moment of this correlator 
was presented in [20], where, aside from the Laporta algorithm, also techniques 
based on Sbases and special treatment of internal self-energies were used for 
checks. These methods will be explained in detail in the paper at hand. Using 
the method mentioned in the previous paragraph, the thirty lowest moments 
of the double-fermionic corrections at 0(a 3 s n 2 ) were determined in (21]. The 

part proportional to 0(a^n\ n ~^) is even known to all orders in perturbation 
theory [22]. 

For the pseudo-scalar correlator, the first moment and second moment at 
0(al) are given in [23], together with the first moment of the axial- vector and 
the scalar correlator. 

In this work, we present the calculation of the second and third moments 
of the vector, axial-vector and scalar correlators and the third and fourth 
moments of the pseudo-scalar correlator. In Section [2] we define our notation 
and explain the methods used in the calculation. In Section [3] we present a 
newly developed method for the efficient treatment of diagrams containing 
self-energy insertions. Section [4] contains the results for the various currents 
and an update of the determination of charm- and bottom-quark masses from 
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experimental data and a s from lattice calculations is given in Section 5. We 
conclude in Section [6j 



2 Methods 

The polarization functions are defined by 

(-?V + »)nV) + q»q»n{( q 2 ) = i J dxe*(o\Tjfa)j' v (o)\o) (1) 

for 5 — v, a, 

q 2 U 5 {q 2 ) = i J dxe iqx (0\Tj 5 {x)j 5 {0)\0) (2) 
for S — s,p, 

with the currents 

jl = $7$, ft = $lnl$ip, f = W>, f = i^l5ip- (3) 
In the low energy limit each of the polarization functions can be written as 

2 

a series in z = ^-y, where m is the mass of the heavy quark and Q q the 
corresponding charge, 

nV) = ^E^» (4) 

n>0 

The coefficients can be expanded in a power series in — : 

C S n = + ^C F C^ + M 2 C^ + (^Y CJW + • • • . (5) 

The decomposition of the four-loop contribution C^' 5 according to the num- 
ber of internal quark loops and its colour structure leads to 

q(3),6 _ CpTpnfC^ 6 + C F T F n 2 h C^ + C F T F n l n h C^ ) J i + C§£ 

+ C F T F rii (C A C^A,n + C F Cf2^ ) + C F T F n h (C A Cj? NAn + C F C^} 5 n ^ . 

(6) 

N 2 -1 

Here C F = 2 c Nc and C A = Nc are the Casimir operators of the fundamental 
and adjoint representation of the SU (N c ) group, respectively Tp = ~ is the 
index of the fundamental representation, rih and n\ denote the number of heavy 
and light quarks, respectively. Co contains the purely bosonic contributions, 

where we set the number of colours Nc = 3 for simplicity. 

As described in [I] the theoretically computed moments C n can be combined 
with measurements of the cross section of e + e~ — > hadrons to determine the 
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masses of charm and bottom quarks. The method is based on the dispersion 
relation 

127r 2 7o s{s-q 2 ) 



n(? 2 ) = T^/„ rf^- (7) 



where 

<r(e + e _ — > hadrons) 
R ( s ) = — rrr 1 — r- 

Taylor expanding both sides of Eq. ([7]) leads to 



with the experimental moments 



Ml xp = / ds 



The calculation of the theoretical moments C n proceeds as follows: The dia- 
grams are generated using qgraf [21]. In total there are 701 four-loop diagrams 
of the propagator type. Some of these diagrams are singlet diagrams, i.e. dia- 
grams with massless cuts, and are not considered in this work. Subsequently 
the diagrams are expanded in the external momentum q 2 and mapped to six 
topologies of vacuum integrals. This is done using q2e and exp [25J in combina- 
tion with Matad [2S] written in FORM [27] . This procedure leads to integrals 
with a maximum of 12 additional powers of propagators and 8 irreducible 
scalar products. The large number of integrals obtained in this step can be 
reduced to a small set of master integrals solving a large system of linear equa- 
tions generated by Integration- by-Parts (IBP) identities [2H]- This reduction is 
achieved using Laporta's algorithm [12J implemented in Crusher [2S]. This 
standard method can be assisted by a special treatment of integrals containing 
self energies, which is described in detail in Section [31 It has to be noted that 
this aforementioned reduction comprises the most difficult part of the calcula- 
tion. The master integrals have been calculated in [30T[3T|l3 2 .33.34 35. 36ll37|l38j . 
After performing the renormalization of the quark masses and the strong cou- 
pling constant in the MS scheme the results given in Section H] are obtained. 



3 Reduction technique for integrals with internal self energies 



3.1 Algorithm 



The required CPU time for the reduction of Feynman integrals of a given 
topology to master integrals using Laporta's algorithm strongly depends on 
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the powers of the propagators and scalar products. For a fixed depth of a Tay- 
lor expansion in an external momentum of a Feynman amplitude the maximal 
number of propagator powers will appear on those integrals which contain the 
maximal number of self energy insertions. In the case of 4-loop tadpoles we 
have at most 3 self energy subgraphs which lead to a sum of propagator powers 
raised by 2 compared to the generic case without self energies. Such integrals 
are thus particularly difficult for Laporta's algorithm. On the other hand the 
presence of self energy insertions can be exploited to perform the reduction 
to master integrals in two less expensive steps, first, the reduction of the self 
energy subgraphs and, second, the reduction of the remaining integral. In the 
second step Intgration-by- Parts relations are constructed in which the self en- 
ergy master integrals are treated as objects depending only on their external 
momentum q. In the following the procedure is illustrated for an arbitrary 
integral which contains a self energy insertion: 




(q,m)PZ(q 2 ) , (11) 



where S / j, 1 ___ IM (q,m) is a rank n tensor self energy integral, R^ 1 - 41n {q 1 m) is the 
rest graph, and Pm {q 2 ) — (q 2 —ml)~ a ° is the connecting propagator with mass 
m . All masses are assumed to be either zero or m. S lil ... ltn (q,m) is expressed 
in terms of a complete set of Lorentz structures {11^ } of rank n consisting 
of the momentum q and the metric g^, and the (scalar) self energy master 
integrals S z (q 2 ,m): 

(q, m) = E 11^ E cP M l\ rn)S z (q 2 , m) . (12) 
p z 

The coefficients c p (d, q 2 , m) are rational functions of the space-time dimension 
d, the momentum q 2 , and the mass m. Inserting this decomposition into ffTTT) 
gives 

T= f rf , ffE^"^(9.^Wi..^(ff«(^)EW,^'n)4(9 2 ) T n ). (13) 

J p z 

The rest graph tensor is contracted with the Lorentz structures to 

R^{q,m)W^ n {q) = R p (q 2 ,m). (14) 

From the structure of the IBP identites we know, that the q 2 dependence of 
the denominators of the coefficients c p z (d, q 2 , m) factorizes into propagator-like 
objects Pk(q 2 ,m) = (q 2 — k 2 m 2 )~ l , where k takes values from a finite set of 
integers depending on the type of the self energy. Partial fractioning leads to 

PZi^yM^rn) = EE^( d ) p r(? 2 >^) ) (15) 

k r k 

5 




Fig. 1. The figure above shows an example of eq. (fT7|) for a four-loop tadpole. Note 
that on the right hand side besides massless (dashed) lines and lines with mass m 
(solid) also a line with mass 2m (double line) appears which is not present in the 
initial integral. The next step would be to repeat the procedure for the two-loop 
self energy. 

where P mo (q 2 ) was written as Pk{q 2 ,m) with m = km. By convention we 
choose r k > for k ^ so that the decomposition (|15p is unique. After 
applying these transformations the integral T takes the form 

T= f A^fffgV)^^^^^)?^? 2 ,™)^^ 2 ,™), (16) 

P z k r k 

which can be written more conveniently as 

T = EEEE %' k ' n (d)T*(ff, r k , m) (17) 

k z V r k 

with 

T z V,r fc ,m) = / d d q R p (q 2 ,m)P r k k (q\m)S z (q\m). (18) 

ff denotes the set of propagator powers of the rest graph. That means the 
initial integral T is expressed as a linear combination of integrals of the type 
T^(ff, rji, m) in which the self energy insertions appear only as master integrals 
and all cross talking momenta between the self energy and the rest graph are 
removed. Figure [Q illustrates eq. (jTSl for a four-loop vacuum diagram. 

The next step is to construct Integration-by-Parts identities for the integrals 
in which the self energy insertions are treated as objects depending only 
on their external momenta. The identities have the form 

= / d d q w £^(q 2 ,m)P^(q 2 ,m)S z (q 2 ,m) (19) 

KP(q 2 ,m)P r k *(q 2 ,m)S z (q 2 ,m) 

I 

= 5 M dT^((f \r k ,m) 

RP(q 2 ,m))p r k *(q 2 ,m)S z (q 2 ,m) 



8 M dT*{f?,r k ,m) + J d d q 







dk, 



+ J d d q 



d 

l dk„ 
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(0,1,2) 



(0,1,3) 



(0,1,2) 



(0,2) 



(0,2) 



(0,2) 



(0,2) 



Fig. 2. All four-loop tadpoles with self energy insertions which appear in our calcu- 
lation can be mapped to these seven topologies. Dotted lines are massless and solid 
lines carry the mass m. The wavy propagator which connects the self energies ap- 
pears with the masses km where k is an element of the list below the corresponding 
diagram. 



where is a loop momentum (i. e. either q or a loop momentum of R p ) 
and i is a loop momentum or, if present, an external momentum. The no- 
tation £ fl [d/dk l j]jR p (q' 2 ,'m) means that the derivative acts on the integrand 
of R p (q 2 ,m). This part and the contribution £^d / 'dk^P^^q 2 ,m) are treated 
like in traditional IBP. The difference is the treatment of d/dk^S z (q 2 ,m) in 
the case k = q where the derivative acts on the self energy. The derivative is 
explicitly performed and the resulting integrals are reduced to the self energy 
master integrals S z . The coefficients appearing in this reduction require partial 
fractioning in q 2 again like in ( 1151) . Afterwards all terms in ( |T9l) are expressed 
by T 2 and the system of equations is solved by the Laporta algorithm. 

In the case of four-loop tadpoles with a one-loop self energy insertion the 
restgraph is a two-loop self energy. Therefore the procedure above is applied 
to the two-loop self energy in a second pass. The topologies which appear in 
our calculation are depicted in Fig. ([2]). 

Let us add a few remarks at this point. The number of integers on which 
a topology T 2 fc depends is given by the number of propagators of the initial 
topology T minus the number of propagators of the self energy subgraphs 
and the cross talking scalar products. In the case of four-loop tadpoles only 
one propagator power of initially ten (eight propagators and two irreducible 
scalar products) survives this procedure. Therefore the combinatorics of the 
IBP system for the T 2 fc is much better behaved than for the initial topology. 
There is a price to pay for this simplification, of course. The reduction of 
subgraphs and the partial fractioning of propagators which carry the same 
momentum but different masses leads to a system of IBP relations which 
couples different topologies. Reducing a single propagator power from a self 
energy can produce two powers on an external leg, therefore in the worst 
case the total number of propagator powers might be doubled. If the rest 
graph is sufficiently simple these issues are easy to deal with. Futhermore, the 
tensor reduction of large powers of cross-talking scalar products will produce 
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huge intermediate expressions and many terms in eq. (ITT]) . These expressions 
simplify to the full reduction to master integrals when the solution of the IBP 
system ( 1191) is inserted. The procedure was implemented in a Mathematica 
program. 

3.2 Implementation in FIRE 

Additionally to the stand-alone implementation mentioned above, for one-loop 
self energy insertions the procedure was also realized with the help of the FIRE 
algorithm [39] . 

Instead of the tensor reduction (j!2p the so-called region-bases feature of FIRE 
is used. The tensor structure arises from scalar products of loop-momenta 
between the self energy and the rest graph. The SBases program [10] can 
construct a basis of recursion relations to reduce the power of these cross- 
talking scalar products to zero without caring about the rest of the intgeral. 
Such a basis is called a region-basis. In the second step the scalar self energy 
integrals are reduced to master integrals. This can either be done by a region- 
basis, or with the help of precalculated reduction tables in combination with 
the Rules feature of FIRE. 

In the case of one-loop self energies with one mass (either massive-massive 
or massive-massless) two master integrals appear, the massive tadpole and 
the self energy with both propagator powers equal to one. Integrals which 
contain the tadpole are treated like usual integrals with one loop less than the 
initial integral, because the tadpole factorizes from the rest of the intgeral. 
The treatment of integrals with the self energy master insertion is a bit more 
complicated. One has to introduce a new object in Feynman diagrams. We 
call this a heavy dot which can be either present (integral has a self energy 
master insertion), or not (integral has a tadpole insertion) and is therefore 
represented by an index which is either zero or one. 




Note that we also have to introduce additional lines which carry masses de- 
pending on the self energy type. The massive- massive self energy needs a 
massless and a double-massive line, the massive-massless self energy needs a 
massless and a usual massive line that is external with respect to the self 
energy (one of these is naturally present in the initial diagram). These lines 
are needed to absorb the momentum dependence of the coefficients from the 
reduction of the self energies to master integrals. The partial fractioning of 
these two lines is done by a recursion relation which is added to the region ba- 
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sis. The IPBs for the remaining integrals are constructed as in eq. (Fl9|) where 
the creation or destruction of a heavy dot is done by a shift operator for the 
corresponding index, just like for usual propagators. The reduction procedure 
aims on reducing the heavy dot index to zero if possible. The FIRE program 
has a setting allowing to use the heavy dots and can create SBases for integrals 
with heavy dots present. However, the IBP generator coming with FIRE at the 
moment is not able to construct such IBP relation. As it has been explained 
above, the procedure can be repeated if a second self energy is present and 
one can again consider a region-basis corresponding to the second self energy 
subdiagram, insert rules reducing it to masters and finally result in a problem 
with two heavy dots. 

3.3 Top level reduction 

Although the approach from above is not applicable to general Feynman inte- 
grals, in combination with FIRE for the top level reduction it can be used for 
the calculation of arbitrary integrals. In the highest sectors the SBases are used 
to reduce the integrals until enough lines are contracted so that only integrals 
with internal self energies are left. From experience we know that the SBases 
algorithm often fails in sectors where self energy insertions are present. At 
this point a program that implements the above procedure is plugged in with 
the rules framework of FIRE to automatically finish the reduction. This way 
we do not need SBases for the lower sectors for which in some cases we could 
not find a basis. In combination with the top level reduction single integrals 
can be calculated if needed or for checks. However, at present this approach is 
not suited for large scale calculations. The self energy formalism itself is very 
efficient and can be used to significantly reduce the effort to invest in Laporta 
algorithm. 



4 Results 

In this Section we present the numerical values of all moments known at 0(a^) 
for the different currents. The full analytical expressions of the new four-loop 
contributions can be found in Appendix |A] We present results up to the third 
moment for the scalar, vector and axial-vector correlators. In the case of the 
pseudo-scalar current the n-th moment is related to the (n — l)-th moment of 
the longitudinal part of the axial-vector current through the Ward identity 

g 2 H£(g 2 ) = 4m\W(q 2 ) - q 2 (dW(q 2 )/dq 2 )\ q i =0 ) . (20) 

This allows to obtain the fourth moment of the pseudo-scalar current. The nu- 
merical values are obtained setting the renormalization scale fi = m and using 
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the MS scheme for the renormalization of the quark masses. We explicitely 
keep the dependence on the number of light quarks n\ and set the number of 
heavy quarks — 1. For completeness we give the numerical values for all 
terms in the perturbative series in a s up to 0(a^). 

cy { cy 

CI =1.06666 + 2.55473— + (0.50988 + 0.66227n,) ( — 

+ (1.87882 - 2.79472n, + 0.09610n z 2 ) {^j , 

CI =0.45714 + 1.10955— + (1.41227 + 0.45491n ; ) — 

7T \1X 

+ (-6.23488 + 0.96156^ - 0.01594nf) (^j , 

(y fa 
Cl =0.27089 + 0.51939— + (0.35222 + 0.42886n,) — 

7T V 7T 

f cy ^ 3 

+ (-8.30971 + 1.94219ra, - 0.03959nf 



7T 



ry f ry 

C\ =0.53333 + 0.84609— + (-2.34665 + 0.41316n,) — 

7T \ 7T 

+ (-1.16280 - 0.56583n, + 0.04784^) (^j , 

ot f d ^ ^ 

C 2 a =0.15238 + 0.14165— + (-0.83002 + 0.19218n,) — 

7T V 7T 



+ (-6.95414 + 1.11092n, - 0.02049nf) 



3 



Cl =0.06772 - 0.01276— + (-0.67592 + 0.13562^) 

7T \ 7T 

+ (-5.36382 + 0.98605n, - 0.02233n? 



a, N 2 



C[ =0.8 + 0.60246— + (-9.50321 + 0.58765n,) (— 



2 



+ (2.36044 - 3.31076^ + 0.23981nf) , 

Oi f Ct 

C s 2 =0.22857 + 0.42582— + (-1.44346 + 0.23664n,) — 

7T V 7T 

fa x 3 

+ (-16.84601 + 1.28345/1; + 0.00398nf) — 

V 7T 

ry / ry \ ^ 

C s 3 =0.10158 + 0.15355— + (-0.60909 + 0.15633n,) — 

7T V 7T 
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+ (-10.40154 + 1.32090^ - 0.01802< 



Ot { OL \ * 

C{ =1.33333 + 3.11111— + (-1.73650 + 0.61728^) — 

TT V TT / 



+ (21.34792 - 8.66336^ + 0.37997^ 



2\ 



a. 



TT 



3 



OL f OL 

Cf =0.53333 + 2.06419— + (6.36704 + 0.28971^) — 

TT V TT 

3 

+ (10.92473 - 1.49687^ + 0.07020nf) ( — 

V TT 

OL f OL 

CI =0.30476 + 1.21171— + (5.19573 + 0.26782^) — 

TT \ TT 

3 

+ (13.96839 + 0.15741^ + 0.01535nf) f — 

Ot f OL 

C\ =0.20317 + 0.71275— + (3.40816 + 0.28627^) — 

TT \ TT 

+ (10.20740 + 1.06763^ - 0.00916nf) (— ^ 

\ TT 



We find full analytical and numerical agreement with all previously known 
results, which means C\ from Refs [T8"1[T9"] and Cf and Cf from Ref. [23J. 

In [6] the value of the third moment of the vector current has been predicted 
C£ = -3.279±0.573, -1.457±0.579 for m = 3,4, respectively. The calculated 
values —2.839 and —1.174 are within these error bound. They also lie well 
within the estimates of [1]. 



5 Applications 



One of the main applications of the low-energy moments of the vacuum polar- 
ization function is the determination of the charm- and bottom-quark masses 
using sum rules. As can be seen from Eq. ([9]) the masses can be extracted 
from the experimental measurement of the i?-ratio in the threshold region 
of open charm or bottom production. This relation has been exploited in [4] 
where only the first moment of the vector correlator was used. Recently, new 
data from the Babar collaboration was published for the region around the bb 
threshold [4T]. Since a full analysis of this new data is beyond the scope of 
this paper, we only quote the results from [42], where besides the new Babar 
data also all moments presented in this work were included. It turns out that 
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moment 


m c (3 GeV)[GeV] Ref. [5] 


m c (3 GeV)[GeVl 


2 


0.986(11) 


0.986(11) 


3 


0.986(10) 


0.986(10) 


4 


0.973(19) 


0.981(13) 


5 


0.969(23) 


0.975(17) 



Table 1 

Update of the results from Table II given in [5] for m c extracted using the pseudo- 
scalar correlator. The value of the fifth moment was taken from Ref. [7]. 

our new results lead to a reduction of the theoretical uncertainty by about 6 
- 10 MeV while the mean values extracted using different moments get only 
shifted by about 2-4 MeV. The values of the quark masses obtained from the 
different moments are in good agreement which demonstrates the consistency 
of the method. The final values for the charm and bottom quark read 

m c (3GeV) = 0.986(13)GeV , (21) 
m b (m b ) = 4.163(16)GeV. (22) 

In the case of the charm quark the experimental data can be replaced by lattice 
"data", i.e. data obtained from the calculation of appropriate correlators on 
the lattice. This has recently been done by the HPQCD collaboration in Ref. 
[5]. They found that in practice the pseudo-scalar correlator is best suited for 
these calculations, but also the other currents can be used with less accuracy. 
In Table [1] we give an update of the results in [5] where our result for the 
third moment of the pseudo-scalar correlator was already used prior to this 
publication. The values of the charm-quark mass extracted using different 
moments are in very good agreement with each other and the results from 
analyses using experimental data as input. 

Lattice calculations can also be used to determine the value of the strong cou- 
pling constant as explained in detail in [5]. Again the pseudo-scalar correlator 
is best suited and in Table [2] we give an update of the numbers presented in 
[5]. The values from different moments are again in very good agreement with 
each other and competitive with the world average for a s [Mz) — 0.1176(20) 



6 Conclusion 

We calculated the second and third low-energy moment of the heavy quark 
correlator for vector, scalar and axial-vector currents and derived the fourth 
moment of the pseudo-scalar correlator. To this end we discussed new meth- 
ods based on Sbases and reduction of internal self-energies for the reduction of 
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moment 


a s (3GeV) Ref. \5\ 


a s (3 GeV) 


a s (M z ) 


1 


0.252(6) 


0.252(6) 


0.1177(12) 


2/3 


0.249(6) 


0.249(6) 


0.1170(13) 


3/4 


0.224(31) 


0.237(11) 


0.1145(25) 


4/5 


0.241(30) 


0.236(19) 


0.1143(44) 



Table 2 

Update of the results for a s from Table II given in [5]. The values of a s are extracted 
using the indicated ratios of moments. The value of the fifth moment was taken from 
Ref. 0. 

scalar integrals to masters. Furthermore we give an update of the results for the 
quark masses and a s obtained in previous publications. The new results pre- 
sented reduce the uncertainty on the quark masses by about 6-10 MeV while 
the central values are only shifted by about 2-4 MeV. This demonstrates 
the validity of the method used for the mass determination. Furthermore the 
new results can be used to improve the reconstruction of n(g 2 ) over the whole 
energy range using Pade approximation. 
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A Analytical Results 

In this Appendix we list the analytical results for the different currents where 
we used the abbreviations 

c 4 =24a 4 + log 4 (2)-6C 2 log 2 (2), 

OO 

a n =Lin(|) = X] ^kt^, ' 
fc=i z K 

OO 1 

Cn = E F - (a.i) 

fc=i k 
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The results are split according to the colour structures defined in Eq. §6§ and 
use the MS-scheme for the renormalization of the quark masses. 

The results of this calculation can be downloaded in computer readable form 
from 

|http : //www-ttp . particle . uni-karlsruhe . d e/Progdata/ttp09/ttp09-18/ 



A.l Second moments 



A. 1.1 Vector current 



c(3)jV _ 15441973 _ 32 

"' 2 19136250 45 

{3)>v _ 1842464707 2744471 

hK2 ~ 646652160 _ 1064448 

{3)}V _95040709 2029 99421 12159109 

lh ' 2 "62705664 ~ 41472 ° 4 + 55296 ^ 4 ~ 4644864 



(3)^ _ 22559166733 520999 167529079 309132631 

lNA > 2 " ~ 16796160000 _ 4354560 ° 4 + 5806080 ^ 4 _ 12902400 ^ 3 ' 
, 3)v 357543003871 520999 598455689 „ 36896356307 , 

lA > 2 11757312000 2177280 4 2903040 S 174182400 ^ 3 ' 
{3)>v 20427854209619 31595849 968787977 362 

hNA ' 2 " ~ 5649153269760 ~ 1161216Q C4 + 15482880 ^ + ~63~^ 5 



q(3),v _ _ : _ = 

hA ' 2 271593907200 2177280 4 1451520 M 

^ 5811074101069 

6706022400 

64985074258811347 1662518706713 26401638588211 



29638030087837 

697426329600 
37320009196157 130387543 2218910663 



_ — ~ ^ ^^-^ 



/n r 2 353072079360000 21016195200 ^ 28021593600 
164928917, 2900811008 1684950406, ..^ 
"C5 ono .,, Qs 0< ., „. c log(2)C 4 



270270 ' 3648645 3648645 
725202752, 362601376, 

log (2)C 2 + log 5 (2) 

10945935 v ; 54729675 v ; 
112680551036302633 

+ 47076277248000 ^ 3 ' 



A. 1.2 Scalar current 



(3)iS _ 381690470169079 66906848^ 402167880157 
nf ' 2 ~ 467026560000 + 1216215 05 + 28021593600 ° 4 
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8363356 , 16726712, 87227759, 

lo g ( 2 ) + orAorAt lQ g ( 2 )C2 - 101P01g log(2)C4 



18243225 ° v 7 3648645 ° v y " 1216215 
1075951565336201 3294508921817 480273679 

560431872000 + 5337446400 + 694980 
(3 ), s _ 355706449933319 4135403 c 857749280491 

hN > 2 " ~ 37661021798400 ~ 2150400° 4 ~~ 37196070912 ^ 3 
83367701 _25 

1720320 21 
( 3)>s _ 102230727187 1121 986734811 4419371 

Z7V ' 2 ~ 19595520000 ~ 48384 ° 4 + 58060800 ^ 3 _ 322560 
(3)jS _ 589278012260869 21022661 98124313871713 

hA ' 2 ~ 2353813862400 + 181440 ° 4 + 58118860800 
357046549 

120960 ^ 4 

(3)ia _ 12729644591 1121 1192980407 5073991 

lA ' 2 ' ~ 279936000 + 24192° 4 + 9676800 ~ 53760 
(3)ja _ 64429199 1301 _ 37067899 63749 

" , ' 2 ~ 41803776 ~ 27648 Q _ 15482880 + 36864 
0),, = 34805101 _ 101467 

hh ' 2 6735960 23760 
p),, = 1439317 _ _8_ 

ll ' 2 6378750 45 



A. 1.3 Axial-vector current 



{3):a _ 217627345572539363 22929814528 48366697030607 

nf > 2 ~ 29957630976000 18243225 ° 5 + 367783416000 ° 4 

2866226816, *, n . 5732453632, 

+ log (2) log 3 (2)C 2 

273648375 v ; 54729675 v ; 

5480723618, , N> 19481983397047414063 _ 

log 2 C 4 Cs 

2606175 6V ;s 3295339407360000 s 

3999530651550337 64302085207 

980755776000 + 12162150 
(3))a _ 58553635489580981 5032917199 c 289668262443149 

fe7V ' 2 " " 3954407288832000 ~ 3048192000° 4 ~ 19527937228800 
40203350213 382 

+ 812851200 ~ 15" 
(3)>a _ _ 414223585133 2041307 43700406067 95269223 

lN ' 2 ' ~ 33592320000 ~ 43545600 ° 4 + 696729600 ^ 3 ~ 1658880 ^ 4 ' 



c (3), a _ 30093199294856077 707672969 c 1650262562573011 

hA > 2 ~ 35307207936000 + 1814400 ° 4 + 290594304000 ^ 3 
1718278879 

172800 ^ 4 ' 

( )a 4855444846081 2041307 30349640587 , 2828510767 , 
lA ' 2 ' 23514624000 21772800 4 49766400 S3 5806080 M ' 



15 



(3), _ 809030351 2797 92264351 137053 

lh ' 2 ~ 627056640 ~ 82944 ° 4 ~~ 46448640 ^ 3 + 110592 
(3))a _ 1173822473 24710027 

M ' 2 ~~ 323326080 ~ 7983360 
( 3),a = 4275743 _ j*2_ 

ll > 2 19136250 135^ 3 ' 



A. 2 Third moments 



A. 2.1 Vector current 



(3)it , _ 31556642272 _ 256 

"• 3 "49228003125 ~ 40l> 

{3)tV _56877138427 6184964549 

hh ' 3 "12609717120 ~ 1556755200 

(3))t , _ 60361465477 _ 1765 86485 _ 57669161 

lh ' 3 "29393280000 ~~ 31104° 4 + 41472 ^ 4 ~ 17418240^ 3 ' 
(3))t , 1475149211788337 8529817 1510937903 

lNA ' 3 ~ ~ 6452412825600000 ~ 774 14400 ° 4 + 14745600 ^ 4 
561258009401 

~ 6193152000 ^ 3 ' 
{3)tV _983812946922223 8529817 21972351293 

lA ' 3 4389396480000 + 38707200 ° 4 + 17203200 ^ 4 
_ 28995540810097 

21676032000 ^ 3 ' 
(3)it , 454880458419083629 7110196837 c 1068488091383 

WVA ' 3 " " 5854170457175040000 _ 11 17670400 ° 4 + 7451136000 
4448 43875740175477222611 

+ ~3ll> 5 ~ 433642256087040000 ^ 3 ' 
(3) lt , _ _ 2327115263308753 16870125343 286864384271 

hA ' 3 " 2489610816000 39916800 ° 4 + 26611200 ^ 4 

377837317054807 

"C3) 



61471872000 

(3))t , _ 801 1001677156303009183663 16091704629458603 

ra °/- 3 1270551061901721600000 45731240755200 ° 4 

1505000915688143609 1781851011826 859399602944 



304874938368000 ' 310134825 ^ 310134825 
38830116184, . x> 214849900736, „ 

" 44304975 l0g(2)C4 " 930404475 log (2)C2 
107424950368, 1061162538194750079871, 

H log(2) 5 H (3 • 

4652022375 &V ; 128047474114560000 S 
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A. 2. 2 Scalar current 



(3)li _ 2229649183 _ _64_ 

ll ' 3 16409334375 405 

(3)ia _ 82738947097 11973270941 

hh ' 3 "12009254400 ~ 2075673600 



iM ~~ 627056640000 " 331 7760 ° 4 + 4423680 ^ 4 _ 1857945600^ 3 ' 
(3)iS 22504447164310009 2832217 288902173 

lNA ' 3 " ~ 1075402137600000 ~ 116 1 2 1 60 04 ~ 4423680 ^ 4 
142693634059 

+ 1857945600 
(3)ia _ _ 296129104756579 2832217 _ 51185293321 

M ' 3 = ~ 1024192512000 + 58060800 04 77414400 ^ 4 
1811694475921 

+ 2167603200 ^' 
(3)ia 549859189388336126083 _ 270366644899 1953161743091 

hNA ' 3 ~ ~ 31222242438266880000 ~ 89413632000 ° 4 + 23843635200 ^ 4 
1888 16019821402882223177 

" ~315~^ 5 ~~ 462551739826176000 
(3)iS _ 1248837940639034917 108697503829 263931412979 

hA ' 3 838406080512000 + 159667200 ° 4 15206400 ^ 4 
4308361964639349413 

+ 434729078784000 ^ 3 ' 
(3)ia _ 1831719443479347906961831 23476362679129111 

"/• 3 211758510316953600000 + 266765571072000 ° 4 

227245780513489691 550981545544 31233089024 

71137485619200 ^ 4 + 103378275 ^ 5 " 103378275 ° 5 
132939852688, , x> 7808272256, 

bg(2)C 4 og 3 (2)C 2 

103378275 SV A 310134825 S V A 

3904136128, 6840706109244803149, 

+ log(2) 5 Cs • 

1550674125 6V ; 798870159360000 s 



A. 2. 3 Axial-vector current 



( 3), a _ 6052378456 64 

"' 3 "49228003125 ~~ 40l> 

(3))a _ 599209514131 _ 12494549257 

hh > 3 "126097171200 ~~ 3113510400 ^ 3 ' 

(3))a _ 4962625369889 630217 30880633 21887209193 

lh ' 3 "3762339840000 " 19906560° 4 + 26542080 ^ 4 " 11147673600 
3))fl 41502247972857197 254785409 2632884019279 

VA ' 3 " " 955913011200000 ~~ 8360755200 ° 4 ~ 11147673600 ^ 4 
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55622383674929 

+ 222953472000 ^ 3 ' 
c (3), a _ 407919130830390283 254785409 17645944577969 

lA ' 3 ' 368709304320000 + 4180377600° 4 5573836800 ^ 4 
2942092239074473 

+ 780337152000 
(3)ia _ 9291278450921361191597 _ 131718995983 

feAM ' 3 ~ ~ 266210698684170240000 ~ 1609445376000 04 

12854292041443 9376 29030831574992389907 

+ 429185433600 ^ 4 ~ lbl> 5 + 1314620734242816000 
( 3), a _ 167314115446749909343 6229030586993 ^ 15127004909623 

hA > 3 35213055381504000 + 2874009600 ° 4 273715200 ^ 4 
49375038643781532767 

+ 1565024683622400 
( 3), a _ 1229252060576583400157229437 26069736016065697819 ^ 

n /- 3 ~ 228699191 14230988800000 + 28810681675776000 ° 4 
84934038473664674083 29147971741346 

1536569689374720 ^ 4 + 930404475 ^ 5 
6290348482816 11153783461592, 

a 5 log(2)C4 

930404475 930404475 v y 

1572587120704, s . _ 786293560352, 

log (2)C 2 H log(2) & 

2791213425 v ; 13956067125 v ; 

_ 18647336379797344185371 
1075598782562304000 ^ 3 ' 



Pseudo-scalar current 



c(3)tP _ 5893303 _ _64_ 

"' 3 9568125 135 

(3)iP _ 26179537 76295 

ftM " ~ 58786560 + 290304^ 3 ' 



103200341 1435 70315 5574817 
627056640 ~ 82944 ° 4 + 110592 ~ 9289728 ( 



(3)i p _ 102523166819 259289 3760967 927288197 

lNA > 3 ~ 11197440000 ~ 4838400 ° 4 + 36864 ^ 4 ~ 9289728 
, 3)rJ 1867646866703 259289 184831769 , 16109807909 , 

lA ' 3 ~ 7838208000 2419200 4 215040 16588800 
( 3)) p _ 10104443767785889 4409031071 25692311767 

AJVA ' 3 ~~ 1318135762944000 ~ 677376000 ° 4 + 180633600 ^ 4 
368 248471102645509 

+ "il> ~ 2169770803200 ^ 3 ' 
(3), P _ 11809333374019073 1657837871 28182137887 

hA > 3 ~ 11769069312000 3628800 ° 4 + 2419200 ^ 4 

214584455705773 

32288256000 
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(3 ),p_ _ 841507212739667387 62193803213161 9141288770183 

n /> 3 ~ 549223234560000 245188944000 ° 4 5108103000 ^ 4 
58401651592 > 15214899712 3649219472, . x> 

C5 H a 5 H loe(2)C4 

6081075 S 6081075 868725 SV A 

3803724928, 1901862464, = 

+ log (2)C 2 log 5 (2) 

18243225 6 v ;s 91216125 6 v ; 

767220341123064149 
+ 84495882240000 ^ 3 ' 



A. 3 Fourth Moments 



A. 3.1 Pseudo- scalar current 



(3)l p _ 329624056 _ _64_ 

"' 4 "607753125 135^ 3 ' 
(3)iP _ 1739021393 52674403 

hM "4670265600 ~ 115315200 ^ 3 ' 

(3)iP _6245190619 18589 910861 1587980983 

lhA "9289728000 ~ 737280 Q + 983040 ^ 4 ~ 1238630400 
(3);P _ 331245076513169857 964223 63056544809 

lNAA 8603217100800000 ~~ 15925248 ° 4 + 148635648 ^ 4 
1229566890599 

-C3, 



2972712960 

(3)iP _ 1309551728160611 964223 10470510368767 

lAA 910393344000 + 7962624° 4 + 1857945600 ( 
36245830807891 

C3, 



5780275200 

(3)) p _ 33124023975801583681957 _ 2883014289719 c 

ftArA ' 4 ~ 1686001091666411520000 ~ 107296358400 ° 4 

87140698714231 11776 3691162910334747972379 

+ 143061811200 ^ 4 + 189 ^ 5 8325931316871168000 ^' 
(3))P _ _ 217480554135087391943 _ 2695980497779 

feA ' 4 ~ 35213055381504000 958003200 ° 4 

45831907340963 106836033862656915617 

+ 638668800 ^ 4 2608374472704000 ^ 3 ' 
(3))P _ 160999779080808137440043699 15601462649447899319 

n /> 4 ~ 7623306371410329600000 9603560558592000 ° 4 

41651006271487731961 3679557604076 133358649856 

12804747411456000 ^ 4 62026965 ^ 5 + 8860995 ° 5 

1559619599504, . ^ 33339662464, „. XA 

+ log 2 C 4 + log (2 C2 

62026965 6V ;s 26582985 6 v ;s 

16669831232, 519180776457371224187 , 

log(2) 5 H (3 • 

132914925 &V ; 21090172207104000 S 
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